Sloshing phenomenon is one of the complex problems in free surface flow phenomena. Numerical meshless methods as a new method can be used to solve this problem. In these methods, the lack of a mesh and complex elements for the domain of problems due to the change in geometry of the solution over time provides a lot of flexibility in solving numerical problems. In the previous researches, the sloshing problem in reservoirs was solved, using the Laplace equation with respect to the velocity potential, but the solution to this problem with pressure equations has not much considered; therefore, using the pressure equations and a suitable lagrangian time algorithm, generalized exponential basis function method has been developed for dynamic stimulation reservoirs. The approximation is solved, using a meshless method of generalized exponential basis functions and the entire domain of problem will discrete to a number of nodes and then with appropriate boundary conditions, the unknowns are approximated. In this study, linear and nonlinear examples have been solved under HARMONIC stimulation, in two-dimensional form of rectangular cube tanks, and the results of them have been compared with the analysis solving methods, other numerical methods, and experimental data. The results show that the present method in two-dimensional mode is very noticeable compared with other available lagrangian methods because of accuracy in solving problem and spending time.